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I. INTRODUCTION 



The goal of this paper is to study possibilities of formulating variational principles for bound- 
ary conditions appearing in the extended thermodynamic systems, where the usual locally con- 
served fields (the mass density, the momentum density, and the energy density) are supplemented 
by various non-conserved fields such as extra stresses in rheology [JT]], higher-order moments of 
the one-particle distribution function in moment systems derived in kinetic theory of gases and 
plasmas, and many others. In order to be specific, we shall restrict our attention to the case of 
so-called extended thermodynamic system, underpinned by Grad's moment method of the Boltz- 
mann equation of rarefied gas Since the seminal work of Grad [@], it is well known that the 
stationary problems in moment equations is ill-posed. Indeed, on physical grounds, it is often un- 
clear how to infer the values of the higher moments on the boundaries without a more microscopic 
considerations. 

In a situation where imposing boundary conditions is problematic or ill-posed, two main di- 
rections in the search for formulations of the boundary conditions can be distinguished. The first 
direction can be broadly characterized as a variational approach. A typical and quite well known 
representative of this strategy are so-called natural variational formulations of stationary equations 
|T3p. This approach is widely used, in particular, in numerical methods based on local minimization 
schemes, such as the finite elements method [Q]. Without going into any detail here, we mention 
that if the solution can be written as a minimizer of a functional, then it is sometimes possible to 
extend the solution from the bulk to the boundary, or to modify the functional in such a way as to 
make this extension possible. By doing so, the natural variational formulations results in so-called 
natural boundary conditions. Many examples are given in the standard references on the finite 
elements method [Q]. It should be also noticed that the physical significance of the boundary con- 
dition thus arising is rarely addressed, especially in the case of extra fields without direct physical 
interpretation. The physics that is behind the behavior in the bulk may not be identical with the 
physics that is behind the boundary conditions. For example, new type of forces arise often on 
boundaries. It is thus possible that a direct extension to boundaries of the potentials that are found 
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to express the physics in the bulk is unrealistic. 

The second strategy is based on an attempt to express the physics that takes place on the 
boundaries. Let us assume that as the result of the physical analysis one formulates a coupled 
system of equations governing the time evolution in both the bulk and on boundaries. States on 
the boundary, or in a boundary layer, are described by values on the boundary (or in boundary 
layers) of the fields chosen to describe states in the bulk and possibly by some other fields defined 
only on boundaries or boundary layers. Some of the boundary state variables are fixed by an 
outside influence. The rest of them, the uncontrollable boundary state variables, evolve in time 
together with the bulk state variables. Let us assume that analysis of the time evolution equations 
shows that the uncontrolled boundary state variables evolve faster than the bulk state variables and 
that they approach, as the time goes to infinity, stationary values. These asymptotically reached 
stationary values of the boundary state variables are then the boundary conditions that we look 
for. We obtained them thus by solving the time evolution equations. If in addition, we are able to 
recognize in the analysis of the fast time evolution a Lyapunov functional, then also this second 
strategy becomes a variational method. This is because the boundary conditions we look for are in 
such a case extremal values of the Lyapunov functional. It is important to emphasize that the way 
the variational functional is introduced in this second strategy does not use the potential arising in 
the bulk, it does not even use the assumption that such potential exists. In fact, it is well known [|5|] 
that the time evolution in the bulk of driven systems can not be often associated with any potential. 
The potential introduced in the second strategy arises from the time evolution of the boundary 
state variables and not from the time evolution in the bulk. 

The second strategy has been mentioned in Ref. [§] as an illustration of a general approach 
to the thermodynamics of driven systems. The potential-driven time evolution of boundary state 
variables have also been used in [[7p in the context of the investigation of consequence of the stick- 
slip boundary conditions in flows of polymeric liquids. The authors of Ref. [|7p do not discuss 
the physical derivation of the boundary time evolution. Also the potential is introduced in [|7|] 
completely phenomenologically. 

Our study in this paper remains also on a phenomenological level. We do not discuss the 
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explicitly the boundary time evolution, we are not therefore in position to recognize the pertinent 
potential in its analysis. We have to use different considerations in order to identify it. Below, we 
shall follow a recent work of Struchtrup and Weiss (see also Ref. [Q]). Struchtrup and Weiss 
[|8p proceed in three steps. 

First, they suggest to consider the local entropy production a as a candidate for the potential 
that will eventually determine the missing boundary conditions. While it is quite well known [|5|] 
that a cannot always be directly related to the time evolution in the bulk, it can still be relevant to 
the boundary conditions (especially in the light of our expectation - based on the physical analysis 
sketched in the previous paragraph - that the boundary time evolution can always be associated 
with a potential). 

Second, having chosen a, one has to ask the question how does this potential depend on the 
boundary conditions. Struchtrup and Weiss [|8|| answer this question as follows: First, they limit 
the analysis to stationary solutions. Let the stationary solution corresponding to a given boundary 
condition is found. The entropy production a, evaluated on the stationary solution, becomes a 
function of both the bulk and the boundary state variables. 

So far, we have arrived at a potential depending on the bulk and the boundary state variables. 
What remains is to make the third step, and eliminate the dependence on the bulk variables. It 
is this third step where our analysis differs from Refs. [||J^]. It has been noticed in Ref. [||] that 
elimination of the bulk variables by averaging the entropy production over the entire volume - and 
which eventually leads to the total entropy production principle in a spirit of Glansdorff and Pri- 
gogine JlOyTTl ] - gives apparently wrong results in applications to the boundary condition problem. 
Instead, a different, much more local analysis has been adopted in [|8],(5p. However, the physical 
significance of such modifications, as well as the physical reasons why the global averaging out 
the bulk variables is not be working have not been addressed. 

In this paper we address the question how physically meaningful variational principles for 
boundary conditions can be constructed on the basis of the entropy production by exploring more 
possibilities than those explored in [0,0]. The intuitive idea behind our consideration is that the 
additional variables in stationary problems often have a significance of a description of the bound- 
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ary layer (this description is greatly reduced, as compared to a full kinetic equation). By adopting 
this viewpoint, we study the question as to what happens if the entropy production is considered 
not in the total volume of the system but rather is localized to sufficiently thin boundary layers. 
A physical interpretation of our results is as follows: If the domain of integration of the entropy 
production is restricted to sufficiently thin boundary layers, the result of the type minimization 
suggests the optimal choice of the boundary condition. Moreover, shrinking the domain where 
the entropy production is sampled from the whole bulk to the boundary layer reveals a behavior 
typical for a critical phenomena, with the optimal value of the boundary condition appearing as 
a result of passing a critical size of the layer. Various features of this transition are studied, and 
plausible realizations of the minimum principle are suggested. 



II. ENTROPY PRODUCTION IN THE BOUNDARY LAYER 

In the context of Grad's method [Q] and its variations, the state of the system is described 
by the locally conserved fields M(x,t) (the local density, momentum and energy), and a finite 
number of nonconserved fields, N(x, t) (nonequilibrium stress tensor, heat flux, fluxes thereof 
etc). The fields N are usually higher order moments of the distribution function which gives a 
full description of the system at a more microscopic level of the kinetic equation. Grad's method 
reduces in a systematic way the description from the level of the kinetic equation for the one-body 
distribution function to the level of a closed set of the moment equations involving only the fields 
M and N. The nonlinear coupled sets of equations in partial derivatives are generically referred 
to as Grad's moment equations, and are given in many sources. The original Grad's method 
[0], technically based on a Hermite polynomial expansion of the one-body distribution function 
satisfying the Boltzmann kinetic equation, has been extended and modified by many authors for 
various kinetic equations [[T^-[T^]. In particular, a generalization of Grad's method to non-moment 



variables has been addressed in [|I3]— |T7|] . Examples of Grad's moment equations will be considered 
in the next section. Here we remind that, each Grad's moment system is equipped with the function 
of the fields, a, the local entropy production. Function cr is nonnegative and equals to zero only at 
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the local equilibrium, and it can be computed once the dissipative terms in the underlying kinetic 
equation are specified (for example, once the Boltzmann collision integral is specified). The form 
of the entropy production also depends on the version of Grad's method used in the derivation 
of moment equations. In many applications, the typical outcome for the entropy production is a 
quadratic form in the fields N (this is valid for small deviations from local equilibrium), 

a = Y,{Ni - (M) ) Aij (M) (Nj - iVf (M)), (1) 

ij 

where iVf^M) are values of the nonconserved fields in the local equilibrium (in terms of the 
kinetic theory, the latter is given by the local equilibrium distribution function which depends 
perimetrically only on the locally conserved fields M, the standard example is the local Maxwell 
distribution function), and where is the positive semidefinite matrix, with matrix elements 
dependent on the functions M and also on the details of particle's interaction in the kinetic picture 
(scattering cross-sections, for example). 

In order to solve the stationary version of Grad's moment equations for the time-independent 
fields, M(x) and N(x), in a domain U E R n , with the boundary dU, a set of boundary conditions 
should be provided. In the typical situation, which we here assume, the boundary conditions for 
the locally conserved fields M (x) are known, and the question concerns only the additional fields 
N(x). To this end, we adopt the first two steps as suggested by Struchtrup and Weiss [§]: First, 
we consider the set of all possible solutions to stationary Grad's equations with the fixed boundary 
conditions for the conserved fields, M b = M(x)\qu, and with various boundary conditions for 
the nonconserved fields, iV b = N(x)\qu- (In principle, other types of boundary conditions could 
be addressed, including derivatives of either M or N, but we shall not consider this option here). 
Second, evaluating the local entropy production functional ([]]) on the configurations of the fields 
thus obtained, we get a set of functions, a(x, N b ), parameterized by the boundary condition data 
N b . Finally, in the third step, one has to eliminate the dependence on x, and to end up with a 
potential ^f(N b ) depending only on the boundary data N b , and whose minima should suggest the 
choice of the boundary condition. It is this third step where we offer a more detailed analysis, as 
compared to Refs. fllU. 
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Specifically, we introduce an additional structure into the domain U . Being inspired by the 
concept of the boundary layer, we introduce a one-parametric family of subdomains B L , where 
L > 0. Each Bl (the test boundary layer) is thought as a subdomain of the characteristic thickness 
L, attached to the boundary dU. For the test boundary layer Bl, we consider the layer-averaged 
entropy production, 

Vol(B L ) Jb l 

where Vol(B L ) is the volume of the subdomain B L . 

The study of minimizers of the set of functionals (Q) for various characteristic thicknesses L 
used to define the boundary layer is the central point of our paper. A priori, it is clear that, if the 
thickness of the layer is taken large enough, then we eventually come to the total bulk-averaged 
entropy production, 

Zu(N b ) = — jf a(x, N b )dx. (3) 

As it has been already demonstrated with explicit examples in Ref. [|8|], minimization of the 
functionals @ over the boundary data N b selects the field configurations beyond a reasonable 
physical interpretation. On the other hand, if we go into the opposite direction, taking thinner test 
boundary layers, and if the hypothesis about the fields N as playing the most important role in the 
description of the physical boundary layer is right, we might expect a crossover in the structure 
of the minimizers of the functional (0). Specifically, we expect that at some value L c , a local 
minimum will start appearing, and which would correspond to the physically plausible value of 
the boundary condition N b . We further expect that variations of these minimal values is not large 
for the entire interval L 6 [0, L c ]. 

This expectation is also motivated in part by the suggestion of Struchtrup and Weiss [|8|] who 
postulated a much more local functional as compared to the total bulk-averaged entropy produc- 
tion (|3|), namely, that correct configurations should minimize the maximum of the local entropy 
production, thus, considering the functional, 

V m (N h ) =maxa(x,N b ). (4) 
xeu 
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Our suggestion to study functionals (Q) which sample the entropy production more locally 
in space as compared to the total entropy production (TJ) does not coincide with the "ultralocal" 
functional (|]), and results are not expected to be identical even in the one-dimensional cases con- 
sidered below. It should be stressed that a correct mathematical definition of the system of the test 
boundary layers B L requires more restrictions but we here do not consider this point rigorously 
here. Finally, the notion of the boundary layer is pertinent to the underlying kinetic theory where 
it can be computed in a few model settings []T8t]. However, it is not straightforward to incorporate 
these results into our considerations. 

Notice that the above construction does not eliminate completely the x-dependence from the 
local entropy production a(x, N b ), rather, it replaces such a dependence by a more transparent one- 
parametric dependence on the thickness L. In principle, any functional T. Bl for L E [0, L c ] can be 
regarded as a potential. However, in practice, a priori estimates for the characteristic value of L 
are sometimes available. These values can be dependent on boundary conditions for controllable 
and uncontrollable fields as well (see next section). On the other hand, the set of the subdomains 
suggests a realization for the potential ^(N b ) which compares the averaged entropy production 
within the layer B L with the averaged entropy production within the rest of the bulk, U \ B L , 

Namely, smoothness of transition from the boundary layer into the bulk suggests the outcome 
for the boundary condition iV b which guarantees that the difference between the averaged entropy 
production in the boundary layer and the averaged entropy production in the bulk is minimal. This 
results in a minimization of the potential, 

^{N b ) = |S BL (iV b ) - Zu\bM)\. (6) 

Other variational principles can be constructed on similar grounds. In particular, if one expects 
that a variation of the entropy production in the boundary layer is considerably higher than in the 
bulk, then the functional (§) can be replaced with a functional involving only local measures 
of activity 
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max a(x, N b 
xeB L 



min a(x, N b ). 

XeB L K ' 



(7) 



Here mm X £B L N b ) represents an approximation to the bulk activity. In the next section we 
shall test all this in the context of a model of heat transfer. 



Following Refs. we consider here a one-dimensional problem of a stationary heat transfer 
for Boltzmann's gas at rest placed between two walls with fixed temperatures. The system is 
described by Grad's 14-moment equations. The set of field variables includes hydrodynamic fields 
M (the mass density p(x), the average velocity v(x), and the temperature T(x)), as well as the 
additional variables N, which are functions of higher moments: The stress tensor t(x), the heat 
flux q(ai), and one more scalar field, A (a?), which corresponds to the fourth-order moment of the 
one-particle distribution function, 



where / cq is the local Maxwellian. 

We further assume that Grad's 14-moment distribution function, f u (M(x), N(x), v), depends 
only on one spatial variable x, and that the velocity dependence is symmetric with respect to 
rotations in the (v y ,v z ) plane. In this case, the average velocity vector, the traceless part of the 
stress tensor r, and the y-, and the ^-components of the heat flux vector, are equal to zero. The 
14-moment Grad's system reduces to the system of four equations for the mass density p(x), for 
the pressure p(x), for the heat flux q(x), and for the fourth moment, A(x), and it reads [|8|J9|JT3|] : 



HI. ONE DIMENSIONAL HEAT CONDUCTION PROBLEM 




d x q = 0, 



d x p = 0, 



d x (A + 15^) = -6 



1 



P 



Up,t) 




i i 



A. 



28 6(P,T) 



(8) 
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Here the pressure p(x) is related to the temperature, and the density as p = (k/m)pT, and the 
positive coefficients £x and £ 2 are the relaxation times, which can be functions of the density p and 
the temperature T. Explicit form of parameters £i j2 is determined by the collision model used in 
the corresponding Boltzmann equation. Following [j^^] we consider two models: the Bhatnagar- 
Gross-Krook equation (BGK model) which gives constant relaxation times £1 = 62 = r, and the 
gas of Maxwell molecules (MM model) which leads to the choice £1 = £2 = 2/3 (ap), where a is 
a constant. The local entropy production for Grad's system (|[) reads: 

_ m 1 J 2 1 1 2 1 1 m A 2 1 

We assume that the walls are placed at x = and x = a. Taking into account the fact that 
the pressure and the heat fluxes are constant, the equations @ require one additional boundary 
condition (in addition to boundary conditions for the temperature T(0) = T , T(a) = T{) in order 
to fix either the heat flux q or the variable A at one of the boundaries. 

Let us first consider the BGK model. It proves convenient to introduce reduced variables, 

T A " 

T = L A' = , q' 



To' p{k/m)T ' " p(%r )W 

' x > aT ° nm 
x = — , cr — — -r — —rO (10) 

Reduced variables ( PD[ ) are used elsewhere below, and we omit primes in order to save notation. 
In terms of variables (|T0|), Grad's equations (|8|) for the BGK model may be written, 

d x (A + 15T) = -G-^—q, 

-''-BGK 

d x T = -- —A, (11) 

28 qK BGK 



where 



■ 1 rr \ V2 
kT \ T 



\m 1 a 

is Knudsen number. The local entropy production @ for the BGK model takes the form, 



1 pg 2 1 A' 2 \ 
a= K^ 5^+120^ • (12) 
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Equations < JTT| ) are easily solved analytically to give 



T(x) = T + W(exp(-x/s) - 1) - 

56 

A(x) = -15W exp(-x/s) H q 2 

5 

Ti - T + 2g/(5ii" BGK ) 




(13) 



exp(— 1/s) — 1 
s = -—qK BGK 



We assume T x > T , then the meaningful values of the heat flux q are negative. Exponential decay 
near the cold boundary x = indicates the boundary layer, and the absolute value of s represents 
its effective thickness (note that if q < then s > 0). Notice that this thickness depends both on 
the Knudsen number, and on the yet unknown boundary condition q. 

As it has been suggested in the previous section, we study minima of the one-parametric family 
of the layer- averaged entropy productions, 



For small Knudsen numbers, and small difference of wall temperatures, results can be com- 
pared with the analytical estimate for q drawn from the conventional Fourier law. In that case, 
as it follows from the Chapman-Enskog solution [|20|], q = — (5/2) K BGK d x T. This allows to an- 
alytically estimate the heat flux as q* ps — (5/2)K bgk (Ti — To). In the test discussed below the 
following set of parameters has been used: T(l) = 1.1, T(0) = 1, and K BGK = 0.05, which results 
in the analytical estimate, q* = —0.0125 for the heat flux, and \s\ w 10 -3 for the characteristic 
thickness of the boundary layer corresponding to this estimate. 

First we compare functionals for various layer widths L. Fig. [I] demonstrates the layer- 
averaged entropy production for different boundary layer thickness L, including the limit 
of the infinitely thin layer, lim^o Si(g), as well as total entropy production, and the functional 
of Struchtrup and Weiss (Q). 

We observe that, when L varies from 1 to 0, there are two qualitatively different outcomes for 
the entropy production E^. For L larger than a crossover value L c , function ( |T4| ) has one 




(14) 
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unphysical minimum q = 0, which coincides with the minimum of the total bulk-averaged entropy 
production (|3|). The latter unphysical minimum has been already reported by Struchtrup and Weiss 
[|8p. However, for L < L c , function E^(g) ( |14"| ) demonstrates another local minimum, q m m(L), 
although the unphysical minimum is still present. As it is seen from Fig. [IJ, variations of the value 
Qmm(L) is small within the interval [0, L c ], and all the values q m i n (L) are close to the analytical 
estimate q*, on the one hand, and on the other hand, these values are close to the minimizer of 
the function (Q). This happens because the maximum of the local entropy production in this and 
similar cases appears to be at the boundary, or within the boundary layer. It is also remarkable that 
there is invariant point where all curves £l(<?) almost touch the curve corresponding to the total 
entropy production. This point is almost the same for any choice of L and it is located very closely 
to the minimum of the function ^i(q) @. 

Fig. |2| compares the three potentials, (@), ^2 (0X an d ^sw @- The value of boundary 
width L in the definition of potentials ^1 and \& 2 was fixed with help of estimate L(q) = s(q) < JT3| ) 
which is the function of boundary condition q. The minima of these functionals correspond to the 
following values of q: 

<hto[*i] = -0.012526 (15) 
g m in[^ 2 ] = -0.012505 
gmin[^sw] = -0.012473 

All these values are very close to the analytical estimate q* = —0.012500. Notice that the estimate 
q* corresponds to the most homogeneous profiles of the local entropy production, and also of 
the temperature (see Fig. [3|, Fig. |], and Fig. |5], |6|). Namely, one observes that if the values q are 
not in small vicinity of q* there is an active domain near left wall x = where an exponential 
decay shows up. It is interesting to note that the boundary layer near right boundary does not have 
any such activity, what is a consequence of the fact that at this boundary the temperature flux is 
directed outward the bulk. In spite of slight deviations in the results obtained with help of different 
potentials they give practically the same temperature profiles. 

Although predictions based on all the three potentials, @, (0), and @, are close to each other 
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in the case of small Knudsen number, we have noticed considerable divergency for larger Knudsen 
number. In order to address this point, we have increased the value of the parameter K BGK , but have 
lowered the value for the dimensionless temperature difference. Figures [7] and || correspond to the 
parameter set K BGK = 0.5, and T\ — T Q = 0.01. We then are able to qualitatively compare this 
result with the direct solutions to the linearized BGK equations reported in the Ref. [|T9p. There 
is a clear indication that when the temperature difference between the walls is sufficiently small 
the solution of BGK kinetic equations gives almost linear temperature profiles in the bulk even 
for large Knudsen numbers. Solution based on our variational principles confirms to this picture 
qualitatively. However Figures 9 and 10 indicate that the Struchtrup-Weiss functional point 
out the solution which is considerably far from "almost linear" unlike the case of small Knudsen 
numbers, which proves that our boundary layer functional are more relevant to the problem of 
selection of boundary conditions. 

Similar analysis has been performed for the model of Maxwell molecules. In terms of variables 
(|T0"1), Grad's moment system for the MM model reads: 

1 q 



d x (A + 15T) = -4 



K T' 



1 1 A 

d x T = -- ——, (16) 

42 K MM qT 

where Knudsen number is, 

(fcTp/m) 3 / 2 

"MM — 

apa 

The local entropy production takes the form, 

1 f q 2 1 A 2 1 

ff =T5xd 4 F+i2T7r <17) 

Because of the nonlinearity, equations ([R]) were solved numerically. For small Knudsen 
numbers, and small difference of the wall temperatures, the heat flux has been estimated as 
q* fa -(15/4)ir MM T'(T 1 - T ), where V = (T + T{)/2. With this, the boundary layer is es- 
timated as L fa ^K MM T'\q\. Like for BGK model we have input the latter estimation into the 
expressions @ and (0) in order to completely specify the functions ^i(q) and ^2(q)- In the test 
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presented below the following parameters were used: K MM = 0.05, T = 1.0, and T\ = 1.1, which 
results in the estimate q* « —0.01969. 

All results for the MM model are similar to those for the BGK model discussed above. Fig. |9| 
demonstrates the crossover in the structure of the layer- averaged entropy production under varia- 
tion of the layer width. Potentials \t> 2 and \l/ sw are compared in Fig. |TD[ Corresponding minima 
of these potentials occur at the following values of the heat flux: 

0.019777, (18) 
0.019714, 
0.019643. 

All these values agree well with the estimate q* w —0.01969. Notice that in both the BGK and 
the MM models, potential ^2 gives the result most close to the analytical prediction. Temperature 
and local entropy production profiles are demonstrated in Fig. [IT] and Fig. [T2[ 

IV. CONCLUSION AND DISCUSSION 

In this paper, we have studied possibilities of introducing a variational principle for boundary 
conditions for Grad moment equations. Our approach is based on a systematic introduction of the 
boundary layer into a phenomenology of variational principles. The approach has been tested for 
models of heat conduction suggested earlier. We have observed that variation of the thickness of 
the domain taken to represent the boundary layer results in a crossover: When L > L c then the 
minimum of the layer- averaged entropy production corresponds to the one predicted be the total 
bulk-averaged entropy production. However, if L < L c , the second local minimum appears, and 
which corresponds to the estimate close to the one resulting from the Struchtrup-Weiss minimax 
principle. This crossover gives an opportunity to define the boundary layer without restoring to 
more precise but also more elaborative microscopic considerations. This observation has led us 
to variational principles which compare the average entropy production in the boundary layer and 
in the bulk. The results have been found in excellent agreement with analytical predictions. The 
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gminf^l] = - 
<?min[^ 2 ] = - 
gmin[^sw] = - 



results of this study therefore make us confident in the usefulness of the entropy production in the 
boundary layer for the problem of boundary conditions in the extended thermodynamic systems. 
The approach is computationally more advantageous than the use of the minimax principle of 
since it avoids a computationally intensive operation of finding extrema of this entropy production 
in entire volume, rather, it is based on a simple integral measure and allow to use simplifications 
for small boundary layer width. 

Finally, it should be stressed that, while the problem of boundary conditions for moment equa- 
tions (and, more broadly, for stationary thermodynamic systems with additional fields) can be ad- 
dressed indeed through consideration of plausible minimum principles, the complete understand- 
ing of those can be accomplished only in the framework of dynamic approach to the boundary 
condition. This point is left for future work. 
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FIGURES 




FIG. 1 . Layer-averaged entropy production £ l (q) ( Jl4| ) as a function of the boundary condition q for 
different layer widths L in the BGK model with K BGK = 0.05 and T\ — To = 0.1. Dashed line is the total 
bulk-averaged entropy production. 
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FIG. 2. Comparison of the potentials (||) and ^ (^) with the Struchtrup-Weiss potential ^sw (0) 
in the BGK model with K BGK = 0.05 and T x - T = 0.1. 
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FIG. 3. Profiles of the reduced temperature T in the BGK model with K BGK = 0.05 and T x - T = 0.1 
corresponding to optimization with various functionals. 
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FIG. 4. Profiles of the reduced local entropy production a (b) in the BGK model with K BGK = 0.05 and 
T\ — Tq = 0.1 corresponding to optimization with various functionals. 
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FIG. 5. The reduced temperature T in BGK model with K B gk = 0.05 and T\ — Tq = 0.1 for large devi- 
ations of boundary condition q from its optimal value q* = —0.0125: curve 1 corresponds to q = —0.005, 
curve 2 to q = —0.01, curve 3 to q = q* (Fourier law), curve 4 to q = —0.015, curve 5 to q = —0.03. 
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FIG. 6. The reduced reduced local entropy production a in BGK model with K BaK = 0.05 and 
T\ — To = 0.1 for large deviations of boundary condition q from its optimal value q* = —0.0125: curve 1 
corresponds to q = —0.005, curve 2 to q = —0.01, curve 3 to q = q* (Fourier law), curve 4 to q = —0.015, 
curve 5 to q = —0.03. 
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FIG. 7. Functionals @ and ^2 (0) as compared to Struchtrup-Weiss potential ^ sw for the BGK 
model with moderate Knudsen number, Kbgk = 0.5 (Ti — Tq = 0.01). 
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FIG. 8. Temperature profiles for moderate Knudsen number, K BaK = 0.5 (Ti — Tq = 0.01), in the BGK 
model. 
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FIG. 9. Layer-averaged entropy production (HH) as a function of the boundary condition q for dif- 

ferent layer widths L in the Maxwell molecules model (If mm = 0.05). Dashed line is the total bulk-averaged 
entropy production. 
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FIG. 10. Comparison of the potential @ and ^ 2 (J7]) with the Struchtrup-Weiss potential ^/ S w (@) 
for the model of Maxwell molecules. 
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FIG. 11. Profiles of the reduced temperature T in the model of Maxwell molecules corresponding to 
minima of various potentials. 
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FIG. 12. Profiles of the reduced local entropy production a in the model of Maxwell molecules corre- 
sponding to minima of various potentials. 
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